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observables associated with the rare decays b — > S7 and b — ► s£ + ^ can be expressed 



In the Standard Model and many of its extensions, it is well known that all of the 
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Q\ ■ in terms of the three Wilson coefficients, C^l^l^qlI/j, ~ nib), together with several 

universal kinematic functions. In particular it has been shown that the numerical values 
^h. of these coefficients can be uniquely extracted by a three parameter global fit to data 



obtainable at future S-factories given sufficient integrated luminosity. In this paper 



we examine if such global fits are also sensitive to new operators beyond those which 
correspond to the above coefficients, i.e., whether is it possible that new operators 
can be of sufficient importance for the three parameter fit to fail and for this to be 
experimentally observable. Using the Left-Right Symmetric Model as an example of a 
scenario with an extended operator basis, we demonstrate via Monte Carlo techniques 
that such a possibility can indeed be realized. In some sense this potential failure of 
the global fit approach can actually be one of its greatest successes in identifying the 
existence of new physics. 
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1 Introduction 



Rare decays of heavy quarks which do not occur at the tree level in the Standard Model (SM) 
can provide a unique opportunity for new physics to reveal itself. When such decays occur 
through loops then the participating SM particles and those associated with the new inter- 
action are placed on the same footing and may yield comparable contributions to the various 
decay amplitudes. In these cases it may be possible to isolate such additional contributions 
and learn something about the detailed nature of the new physics scenario. 

Amongst the rare decays involving b quarks, two of the cleanest and most well under- 
stood inclusive processes are b — > s'j and b — > s£ + £~ . At present, the branching fraction for 
the b -> S7 mode has been measured by CLEO0 to be B(b -> sy) = (2.32±0.57±0.35) x 1(T 4 
while a preliminary result from ALEPHjJ yields the value (3.38±0.74±0.85) x 10~ 4 . On the 
other hand, there exist only upper bounds for the decay b — > s£ + £~; the strongest constraint 
at present is the 90% CL limit B(b -> s£ + £-) < 4.2 x 10- 5 @ from CLEO, which is obtained 
by combining both their di-electron and di-muon data samples. As we will see below, this is 
only a factor of ~ 6 above the expectations of the SM for this branching fraction so that we 
may expect this decay to be observed in the near future. 

In the SM and in many of its extensions (including, e.g., fourth generation mod- 
els, models with an extra down-type quark, SUSY, extended Higgs sectors, Z' scenarios, 
theories with large anomalous gauge boson couplings, etc.) the phenomenology of both of 
these rare processes above are almost completely determined by the numerical values of the 
Wilson coefficients of only a small set of operators evaluated at the scale \x ~ m&. In our 
somewhat unconventional notation these are denoted as CVl,9l,iol(/-0- At the weak scale 
the operators corresponding to these coefficients arise in the SM from the usual 7 and Z 
penguins as well as W box diagrams. It has been successfully argued in the literature M 
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that by combining observables associated with both the b — > 37 and b — > s£ + £~ processes, a 
model independent three dimensional global fit can be performed to numerically determine 
the values of these three Wilson coefficients. Indeed, given sufficient statistics at future 
5-factories this approach leads to only rather modest uncertainties in the fitted values of 
these coefficients allowing us to test the SM and look for new physics. We note that only 
the magnitude |CVl(/x)| can be extracted from the b — > 57 transition so that its sign would 
remain undetermined from this channel alone even if infinite precision were available. For 
the observables associated with the b —> s£ + £~ decay all three of the coefficients contribute 
and therefore their relative signs as well as their magnitudes can be extracted from the data 
when combined with our knowledge of B(b — > sj). 

In some ways the determination of the these three Wilson coefficients via a global 
fitting procedure in rare B decays is similar to the searches for new physics in precision 
electroweak measurements [|] through the use of the oblique parameters S, T, {70. As the 
reader may recall, in the SM (for a reference value of the top quark and Higgs boson masses) 
these parameters are all identically zero. For certain classes of new physics, such as a 
fourth generation of quarks and leptons, fits to precision data would then lead to some 
consistent set of non-zero values for these parameters with a good x 2 . Of course, if these 
parameters are found to be non-zero and it is also found that different precision observables 
yielded statistically distinct S, T, U values then we would necessarily conclude that the new 
interactions are not describable by the oblique corrections alone. (As is well known, any 
class of new interactions that induce significant flavor dependent vertex corrections, such as 
a Z', can not be portrayed solely in terms of S,T,U.) Such a situation would provide a 
unique window on the complex nature of the new physics scenario. 

It is then obvious that when sufficient statistics become available in the not too 
distant future for this type of analysis of rare B decays to be performed, there are only three 
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possible outcomes: (i) The numerical values for the coefficients are found to agree with the 
SM expectations with a good x 2 ■ In this case the new physics is decoupled and either higher 
precision data or searches elsewhere are necessary to uncover it. (ii) A quality fit is obtained 
but the values one finds for the three Wilson coefficients are far from the SM expectations 
in x 2 - This is the result usually discussed and anticipated in the literature f§] in the set 
of extended models listed above, (m) As with the case of precision measurements, the last 
possibility is potentially the most interesting and the one we are interested in here: the value 
of x 2 f° r the best three parameter fit is found to be very large and cannot be accounted for 
by an under estimation of systematic uncertainties. This represents in some sense a failure 
of the model independent approach in that it is clear that the true numerical values of the 
three Wilson coefficients are not being extracted from the data. However, another point of 
view is that this result is in fact this approach's greatest triumph since it is telling us that 
the new interactions necessarily involve an extension of the operator basis to include new 
operators beyond the usual set. This implies that the new physics scenario is richer than 
any of those in the list above. 

Thus the question we wish to address here is whether new physics which does involve 
new operators in an extended basis can indeed manifest itself as a poor fit when we have 
the freedom to vary the 3 coefficients to obtain a good x 2 ■ (Of course, |CVz,(/i)| cannot be 
freely varied by too large an amount due to the reasonable agreement between the present 
data and SM expectations for the b —>■ 57 decay rate as discussed below.) The purpose of 
this paper is then to demonstrate this result by providing an existence proof that a new 
physics scenario of the desired type not only exists, in the form of the Left-Right Symmetric 
Model (LRM), but that it leads to poor values of x 2 m the global fit when only the usual 
three operators are employed. The point we wish to stress here is not the particular physics 
of the LRM, or any other specific model, but that the existence of an extended operator 
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basis can indeed manifest itself in the poorness of the three parameter fit given reasonable 
integrated luminosities. We note, however, that without further analysis the failure of the 
fit itself will not yield information on which new operators would need to be introduced. We 
further note that the LRM is of course not the only new physics scenario with an extended 
operator basis |7]]. 

The outline of this paper is as follows. In Section 2 we overview the status of the 
various pieces necessary for calculations of the b — > 37 and b — > s£ + £~ decay rates and 
distributions in the SM. In Section 3 we provide a background on the basics of the LRM 
and the parameters it contains which are relevant for the processes of interest here. We 
discuss and generate several possible forms of the right-handed CKM weak mixing matrix, 
Vr, following a systematic approach that maintains unitarity and wherein agreement with 
experimental constraints can be easily accommodated. Next we set up the LRM calculations 
for these decays and demonstrate that many LRM parameter space regions exist wherein the 
rate for b — > sj is essentially the same as in the SM. The general formulae for the b — > s£ + £~ 
double differential distributions are discussed and it is shown that the LRM yields distinct 
results for observables associated with this decay even when the b — > sj rate duplicates 
the SM expectation. In Section 4 we discuss our Monte Carlo approach and generate a 
large number of data samples corresponding to each of the several models discussed in the 
previous section. We demonstrate that for high luminosities, corresponding to 5 x 10 s BB 
pairs, typical of samples available at hadron colliders, the resulting fits to the conventional 
three Wilson coefficients lead to very large \ 2 values which clearly signal the failure of the 
fit. For lower luminosities, now corresponding to 5 x 10 7 BB pairs, typical of samples to be 
available at T(4S) machines, the x 2 values are also found to be quite large in most, but not 
all, cases. A discussion of these results and our conclusions can be found in Section 5. 
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2 Rare Decays in the Standard Model 



In order to be convinced that new physics has indeed been discovered it is necessary that the 
predictions for the rates and other observables associated with these decays in the SM be 
on firm ground. In particular, calculations of the decay rates and, in the case of b — ► s£ + £r , 
kinematic distributions have become increasingly sophisticated within the SM context. A 
"straightforward" next-to-leading order(NLO) calculation finds B(b — > sj) = (3.28 ±0.33) x 
10~ 4 [|J. However the inclusion of l/m^ and 1/m^ corrections || increases this value by about 
3%. A further enhancement of about 3% occurs when one systematically disregards any of the 
NNLO terms |TIJ. While closer to the central value of the preliminary ALEPH measurement 
these predictions are certainly consistent with the present CLEO data. For purposes of 
simplicity, in our numerical analysis below we make direct use the NLO result and ignore 
these additional small correction terms. This approximation will have no impact on our 
conclusions. One may anticipate that in the next few years this theoretical uncertainty may 
shrink to as low as 5% as the various input parameters are better determined. A comparable 
experimental determination of this branching fraction may also be possible at future B- 
factories since the measurements will be limited only by systematics. We will assume below 
that the SM value for this branching fraction is essentially realized by future experiments 
within the experimental and theoretical uncertainties. 

In the case of b — > s£ + £~, a complete short-distance NLO calculation has been avail- 
able for some time [11]] and the l/mf correction terms are also known [12| but provide only 



very small modifications. One of the remaining difficulties is the inclusion of non-perturbative 
long distance pieces associated with the J /if} and if)' resonances and the corresponding 
corrections. Here some modeling uncertainties remain and the traditional approach has 
been to treat the resonance contributions phenomenologically, which is not without some 
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difficulties! HI- However, at least in the regions sufficiently below and above the resonances 
(i.e., s = q 2 /ml < 0.3 or > 0.6, where q 2 is the invariant mass of the lepton pair) Buchalla, 
Isidori and Rey|| have shown that the heavy quark expansion in terms of leads to 

reliable predictions with only small corrections to the NLO results and that there are no 
difficulties associated with double counting. In the same kinematic regions the phenomeno- 
logical resonance models give comparable numerical predictions. In many extensions to the 
SM, the l/m 2 b and resonant contributions are either of the same form as in the SM or can 
be easily obtained via suitable modifications of the SM terms and we employ these results 
in our analysis below. We note that the SM predicts a branching fraction for b — > sfi + fi~ of 
~ 6 x 10~ 6 , which is not too far from the present upper bound. 

In order to obtain the complete parton level NLO predictions for these two processes 
in the SM (or in other models), several steps are necessary. First, the complete operator 
basis must be determined at the high (matching) scale, typically taken to be Mw- Secondly, 
the matching conditions for the coefficients of the operators at the high scale must be calcu- 
lated at both the LO and NLO. Thirdly, the anomalous dimension matrices for the relevant 
operators at both LO and NLO are determined and the coefficients are evolved to the fi ~ mj 
scale via the Renormalization Group Equations. Lastly, at the scale fi the matrix elements 
of the relevant operators need to be computed through NLO. For the SM all of these pieces 
are now essentially in place for both the b —>■ s'j and b — > s£ + £~ decays after an enormous 
amount of labor. Unfortunately, the corresponding results only partially exist for most of 
these pieces in almost all extensions to the SM[|14j|. 



3 The Left-Right Model 
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3.1 Model Background 



In order to be self-contained we briefly review the relevant parts of the LRM needed for the 



discussion below; for details of the model the reader is referred to Ref . |15| . The LRM is 
based on the extended gauge group SU(2) L x SU(2) R x U{1) and can lead to interesting 
new effects in the B system |TE[]. Due to the extended gauge structure there are both new 
neutral and charged gauge bosons, Z' and W R , in addition to those present in the Standard 
Model. In this scenario the left- (right-) handed (LH, RH) fermions of the SM are assigned to 
doublets under the SU{2)l( R ) group and a RH neutrino is introduced. The Higgs fields which 
can directly generate SM fermion masses are thus in bi-doublet representations, i.e., they 
transform as doublets under both SU{2) groups. The LRM is quite robust and possesses 
a large number of free parameters which play an interdependent role in the calculation of 
observables and in obtaining the existing constraints on the model resulting from various 
experiments. 

As far as B physics and the subsequent discussion are concerned there are several 
parameters of direct interest, most of which result from the structure and spontaneous sym- 
metry breaking of the extended gauge sector. The ratio of the SU{2) R and SU{2) L gauge 
couplings is bounded by 0.55 < k = g R /gL < 2, where the lower limit is a model constraint 
and the upper one is simply a naturalness assumption. Whereas g^ is directly related to e as 
usual through sin 2 ^, g R is unconstrained except through the definition of electric charge 
and naturalness arguments; GUT embedding scenarios generally suggest that k < 1. For 
simplicity we assume k — 1 in most of our discussion below. The SU{2) L x SU{2) R x U{1) 
extended symmetry is broken down to the SM via the action of Higgs fields that transform 
either as doublets as discussed above, or also possibly as triplets under SU{2) R . This choice 
of Higgs representation determines both the mass relationship between the Z' and W R (anal- 
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ogous to the condition that p = 1 in the SM) as well as the nature of neutrino masses. In 
particular, the Higgs triplet choice allows for the implementation of the see-saw mechanism 
and yields a heavy RH neutrino. We assume triplet breaking below so that the Z' mass is 
completely specified by the Wr mass and the value of k. 

After complete symmetry breaking the Wl and Wr bosons mix; this mixing being 
described by two parameters, a real mixing angle <ft an d a phase u. Note that it is usually 
ttj, = tan (ft which appears in expressions directly related to observables. The additional phase, 
as always, can be a new source of CP violation. However, in discussing processes in which 
the RH neutrinos do not participate, as is the case in B decays, this angle can be thought of 
as an overall phase of the right-handed CKM matrix, Vr, and we will subsequently ignore it. 
The mixing between Wl and Wr results in the mass eigenstates W\ t 2, with a ratio of masses 
r = M\jM\ (with M 2 — Mr). In most models is then naturally of order a few times r, 
or less, in the large M 2 limit. Of course, W\ is the state directly being produced at both the 
Tevatron and LEPII and is identical to the SM W in the limit <fi — > 0. We note that when <fi 
is non-zero, W\ no longer couples to a purely LH current. Of course if a heavy RH neutrino 
is indeed realized then the effective leptonic current coupling to W\ remains purely LH as 
far as low energy experiments are concerned. As is well-known, one of the strongest classical 
constraints on this model arises from polarized /i decay [17], which is trivial to satisfy in this 



case. 



It is important to recall that the extended Higgs sector associated with both the 
breaking of the full LRM gauge group down to U{l) em and with the complete generation of 
fermion masses may also play an important role in low energy physics through the existence 
of complex Yukawa and/or flavor-changing neutral current type couplings. However, this 
sector of the LRM is highly model dependent and is of course quite sensitive to the detailed 
nature of the fermion mass generation problem. For purposes of brevity and simplicity 
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these too will be ignored in the following discussion and we will focus solely on the effects 
associated with W\p exchange. 

Additional parameters arise in the quark sector. In principle the effective mass ma- 
trices for the SM fermions may be non-hermitian implying that the two matrices involved in 
the bi-unitary transformation needed to diagonalize them will be unrelated. This means that 
the elements of the mixing matrix Vr appearing in the RH charged current for quarks will be 
unrelated to the corresponding elements of Vl = Vckm- Vr will then involve 3 new angles 
as well as 6 additional phases all of which are a priori unknown parameters. Needless to say 
the additional phases can be a further source of CP violation. The possibility that Vj, and 
Vr may be unrelated is often overlooked when considering the potential impact of the LRM 
on low energy physics and there has been very little detailed exploration of this more general 
situation. Clearly as the elements of Vr are allowed to vary the impact of the extended gauge 
sector on B physics will be greatly affected. Other well-known constraints on the LRM, such 
as universality, the apparent observed unitarity of the CKM matrix, B° — B° mixing, the 
K L — K s mass difference ^ ], as well as Tevatron direct W searches are quite sensitive 



to variations in Vr[[2(|, but W 2 masses as low as 450 — 500 GeV can be easily accommodated 
by the present data. To be safe and to keep future W 2 searches in mind, however, we will 
generally assume that M 2 > 800 GeV for any form of Vr, implying that the magnitude of 
is less than a few ■ 10~ 2 . In our numerical study below we will assume various different 
forms for Vr; in all cases we will assume for simplicity that the values of the elements of Vl 
as extracted by current experiment |2T| are not much influenced by the existence of the new 
LRM interactions. An updated analysis on the possible general structure of Vr has yet to 
be performed. 
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3.2 Forms of V R 



In order to study the potential influence of the LRM on B physics various forms of Vr 
should be examined. In fact, the possibility that Vl 7^ Vr and that B decays are purely 
right-handed was entertained some time ago by Gronau and Wakaizumi||22||. Just how large 



the right-handed contribution to b — > c transitions is allowed to be within the LRM context 



is not yet accurately known p3| but may be sizeable in magnitude with an unknown relative 



phase. The experimental bounds are as follows. L3 has compared^] their measurements 



of both the lepton and missing energy spectra in semileptonic b — > c decays with a number 
of different hypotheses, and have excluded both the (V + A) x (V — A) and V x (V — A) 
scenarios, clearly indicating that this coupling is dominantly left-handed. This qualitative 



result is confirmed by the sign of the A& polarization observed by ALEPH[^5| at the Z pole. 
The strongest constraint comes from CLEO|26|], with measurements of both the leptonic 
forward-backward asymmetry as well as the D* polarization in the decay B — > D*£u. On 



the theoretical side, Voloshin|]T6[ has recently considered how significant right-handed b — > c 
currents, at the ~ 15% level, may assist in our understanding of the B semileptonic branching 
fraction. In a more general context, other forms of Vr have been discussed by Langacker 



and Sankar|27] upon which we generalize in the following analysis. 



Since we are not concerned with CP violation in what follows, for numerical simplicity 
we will ignore the phases in this matrix in which case it can be completely described by three 
mixing angles. Even in this limiting scenario the set of possible forms for Vr is enormous; 
however, it is sufficient for our purposes here to simply examine some sample forms. We will 
assume that each row and column of Vr contains only one large element with a magnitude 
near unity as is the true for the conventional CKM matrix. In this case there are only six 
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matrices about which we can perturb; we write these symbolically as 



M, 



1 
1 
1 



M 



D 



1 

1 

1 



M 



13 



1 

1 
1 



M, 



1 

1 
1 



M 



c 



1 
1 
1 



1 

1 

1 



For example, matrix D corresponds, in the standard CKM parameterization ||28|| , to the 



situation where si 2 , s 2 3 and Ci 3 are all ~ 1. Following Wolfenstein|29|], this suggests taking 
C12 ~ A n , S13 ~ AX m , and C23 ~ BX P , where A, B are order unity, n,m,p> 1 and A ~ 0.22. 
Of course, all these parameters are not arbitrary since the experimental constraints from 
Kl — K$ and B — B mixing must be satisfied. For r < 10~ 2 and < a few x 10~ 2 , these 
bounds are trivially fulfilled, without further fine tuning, if p > 3 and n, m > 1. For p > 3, 
B physics is not very sensitive to the value of m so we take for simplicity p = 3 and m = 1. 
This leaves us with a set of matrices we can label as Vr = D(n), for n=l,2,3,.... Table 1 
lists the complete set of parameterizations for all six types of Vr mixing matrices, which we 
arrive at by following a similar procedure, and the constrained range of the exponents we use 
in the analysis below. The values of these powers reflect simplicity as well as that required 
to satisfy the low energy experimental constraints. 

Note that we have not employed any constraint on the allowed strength of the right- 
handed b — >• c coupling in our discussion. Since these bounds are relatively complicated their 
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Matrix 


Sl2 


Cl2 


Sl3 


C13 


S23 


C23 




constraints 


A 


X n 


1 


BX P 


1 


AX m 


1 


n 


= 2,p = 3,m = 1,2,3 


B 


1 


X n 


B\ p 


1 


AX m 


1 


n 


= 2, m = 3,p = 1, 2, 3 


C 


AX m 


1 


BX P 


1 


1 


A n 


P 


= l,m = 2, n = 1,2,3 


D 


1 


A n 


AX m 


1 


1 


BX p 


P 


= 3, to = 1, n = 1, 2, 3 


E 


AX m 


1 


1 


A" 


1 


BX P 




p, m > 2, n > 1 


F 


AX m 


1 


1 


A n 


BX P 


1 




p, TO > 2, 72 > 4 



Table 1: Parameterizations of Vr in the absence of CP violation assuming A ~ 0.2, A and 
B are of order unity, and m,n,p > 1. 



detailed numerical impact will be discussed elsewhere [f23|], but will have the most impact in 
the cases of matrices C and D. 

3.3 Rare Decay Formalism 

The analysis of the decays b — > 57 and b — > s£ + £~ in the LRM begins with the following 
extended effective Hamiltonian, 

AC 12 

^// = E CMO iL {^) + L^R, (2) 

The OtL R are the complete set of operators involving only the light fields which govern b — > s 
transitions. The first thing to notice is that for convenience the conventional CKM factors 
in front have been absorbed into the values of the coefficients themselves. This is generally a 
useful approach to follow when multiple mixing angle structures appear simultaneously. The 
second point to note here is that whereas there are only 10 local operators describing b — > s 
transitions in the SM, i.e., On — Oiql, here there are 24 operators. 10 of the additional 
operators correspond to the chiral partners, with L — > R, of those present in the SM . The 
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complete basis for each helicity structure then consists of the usual six 4-quark operators 
@i-6L,r, the penguin-induced electro- and chromo-magnetic operators respectively denoted 
as 7:8LjR , as well as O gLiR ~ es LjRa ^b L>Ra £jH, and O WLjR ~ esL^^nJ^d which 
arise from box diagrams and electroweak penguins. (Here the indices a label the color 
structure of the operators.) In the LRM we not only have the augmentation of the operator 
basis via the obvious doubling of L — > R, but two new additional 4-quark operators of each 
helicity structure are also present at the tree-level due to a possible mixing between the 
W l ,r gauge bosons: 1X l,r ~ {sal^R^cp^b^^R and 12L ,r ~ {sal^R^c^bp)^. 
Operators Oul,r occur at the tree level in a fashion analogous to operators 2 l,r- 

In evolving down from the weak scale to fi ~ rrib these operators mix under renor- 
malization as usual. Fortunately the 24 x 24 anomalous dimension matrices split into two 
identical 12 x 12 chiral submatrices since the operators of each chirality do not mix under 
RGE evolution. The complete 12 x 12 anomalous dimension matrix at LO was first calcu- 
lated by Cho and Misiak|3(J and at NLO only the 10 x 10 submatrix corresponding to the 
SM operators is presently known. 

The determination of the matching conditions for the 24 operators at the electroweak 
scale even at LO is already somewhat cumbersome since the LRM contains a large number of 
parameters and, in addition to new tree graphs, 116 one- loop graphs must also be calculated. 
(Additional diagrams due to possible physical Higgs exchange are not yet included and 
would introduce additional model dependence.) Some of these diagrams have already been 
calculated for the earlier analyses of the decay b —>■ 57 in the LRMf30|j. Clearly, NLO 



matching conditions do not yet exist for this model so we employ only the LO ones below. 
However we note that the numerical size of these NLO contributions in the SM case were 
found to be small. In addition one has to separately include the new LRM contributions 

to the semileptonic branching fraction, Bi, including both finite m c jm^ ~ 0.29 and LO 
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QCD corrections, since it is conventionally used to normalize both the b — > 57 and b — > 
s£ + £~ decay rates. (We assume that the relevant elements of have the same numerical 
values in the LRM as in the SM in our numerical calculations; this need not be the case 
experimentally [p3[ | .) In the case of b — > 57 the required NLO real and virtual corrections 
to the operator matrix elements for the LRM are almost completely obtainable from the 
SM results when augmented by the new terms in the LO anomalous dimension matrix and 
through the use of left -aright symmetry Additional terms arising from operators 0\2l,r have 
yet to be included. Clearly we cannot claim to be performing a complete NLO treatment of 
the LRM case until all of the missing pieces have been calculated. However, this complete 
calculation is not necessary to demonstrate our points but certainly our analysis should be 
repeated once all the NLO pieces are in place for the LRM. The advantage of the present 
approach, however, is that when we turn off the effects associated with the various LRM 
contributions to both b — > 37 and b —>■ s£ + £~ we recover the usual NLO SM results. 

3.4 b -> S7 in the LRM 

Using the central values of the quantities m t (m t ) = 167 GeV, \i = nib, ol s {Mz) = 0.118, 
m c /md = 0.29, k = 1 and Bi = 0.1023 we can calculate the rate for b — > 57 in the LRM 
using the existing pieces of the NLO calculation that are presently available if we also provide 
sample values for the quantities Mw 2 , and the relevant elements of the matrix Vr. Recall 
that we are looking here for particular non-decoupling regions in the LRM parameter space 
that essentially give the same result as the SM for the b — > 37 branching fraction. (We are 
not interested in parameter space regions which differ only infinitesimally from the SM.) Let 
us first examine the simple scenario where Vl = Vr\ this is the so-called manifest LRM. 
In this discussed above, the — K$ mass difference and direct Tevatron collider 

searches require that W 2 be heavy; for purposes of demonstration we take M W2 = 1.6 TeV 
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so that ttj, = tan0 is now the only free parameter as the W2 contributions are now almost 
completely decoupled. Here we note that one of the interesting features uncovered by earlier 



analyses of b — > 57 in the LRMpOf was that left-right mixing terms associated with ^ 
can be enhanced by a helicity flip factor of ~ mt/rrib and can lead to significantly different 
predictions than the SM even in this case where Vl = Vr and W2 is very heavy. 
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Figure 1: Prediction for the b — > S7 branching fraction for the case Vl = Vr and M^ 2 = 1.6 
TeV as a function of t in the LRM. The 95% CL CLEO(ALEPH) allowed range lies inside 
the dashed(dotted) lines. 

Figure [l] shows the prediction for the b — > 57 branching fraction in this case and we 
see that the SM result is reproduced when = 0, as expected, apart from a very small 
correction of order M^ /M^r . However we also see that a non-decoupling conspiratorial 
solution occurs when ~ —0.02, which yields a result which is exactly the same as the SM. 

15 




From using b — » s'y alone, the SM and this LRM case are indistinguishable, independent 
of what further improvements can be made in the branching fraction determination. This 
means that additional measurements would be necessary to separate these two cases. Of 
course, given present experimental data both the ALEPH and CLEO results allow for a 
rather broad range of t^. 

We can now ask if corresponding conspiratorial regions exist for any of the more 
general forms of Vr considered above. For example, if we calculate the b — > 37 rate with the 
Aim) and D{n) matrices we obtain the results in Fig. ^ For the D{n) case, we see that only 
n = 2 provides such a region. For the other matrices we find conspiratorial regions when 
m = 1,2,3 for matrix A and p = 2, 3 for matrix B; no such regions are found for matrices 
C, E and F. Note that m = 2 in case A essentially corresponds to Vl = Vr. We denote 
these six conspiratorial cases as Vl = Vr, A(l,3), 5(2,3) and D(2) for later purposes. It is 
interesting to note that all of these cases lead to modest increases |3T| in B(b — > sg) by as 



much as 75% above the SM prediction at LO. Of course, in addition to these conspiratorial 
solutions, we note from Fig. ^ that a rather wide range of is again allowed by current 
data from CLEO and ALEPH. 

3.5 b -> s£ + £- in the LRM 

For b — > s£ + £~ , the effective Hamiltonian leads to the matrix element (neglecting the strange 
quark mass but keeping the mass of the leptons) 



, , \P1Gfol 
M = 



7T 



(3) 



-2C e 7 f L f m b s L ia^b R hn + L^R 
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Figure 2: Prediction for the b — > 57 branching fraction for the cases Vr = A{m) (top) and 
V R = D(n) (bottom) with k — 1 and M W2 =0.8 TeV as a functions of ^ in the LRM. The 
solid, dashed, dash-dotted curves correspond to m,n=l, 2 and 3, respectively. As before, 
the 95% CL CLEO(ALEPH) allowed range lies inside horizontal the dashed(dotted) lines. 
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where q 2 is the momentum transferred to the lepton pair. Note that Cg{f R contains the usual 
phenomenological long distance terms and that all the CKM elements are now contained in 
the coefficients themselves. From here we can directly obtain the expression for the double 
differential decay distribution 



dB 

dzds 



BiK^-zBtl - s) 



[a 2 _ + a\) + (b 2 _ + b\ 



1 + s) - (1 - s)f3 2 z 2 



(at - a + ) - (b_ - b + ) (3zs + 4x(a_a + + 



+-AC' 7L + C 2 R )(l-s)\l-p 2 z 



(4) 



-Re [C 7L (a_ + a+) + C 7R (b^ + b + )] (1 - s)(l - f3 2 z 2 ' 
s 



where z = cos^+^-, s = q 2 /m 2 , x = mj/m 2 , [3 = Jl — 4x/s, a± = ±Ci l + 2C 7 l/s and 

b± = a± with the replacement L —>■ R. As usual 6e+e- is the angle between the lepton direc- 
tion and that of the original b quark in the lepton pair center of mass frame. Multiplication 
of the above expression by the total number of produced B mesons gives the double differen- 
tial event distribution. We normalize this rate to the usual semileptonic branching fraction 
(Bi = 0.1023), including finite m c jm^ = 0.29 and QCD corrections with a s (Mz) = 0.118, 
which are fully included in the overall normalization parameter K (see Ali et al, in Ref. 0]). 
LRM corrections to the semileptonic rate are, of course, also included; here the assumption 
that the mass of the right-handed neutrino tjin > rrib becomes relevant. For numerical 
purposes we take = 250 GeV in the calculation of the right-handed box diagram for 
b — > s£ + £~, but the precise value oim^ is not important for our purposes of demonstration. 



From this double differential distribution we can compute both the lepton pair in- 
variant mass distribution by integration over z as well as the leptonic forward-backward 
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asymmetry. The asymmetry is given by the expression 

These two observables are shown in Fig. ^ for the SM as well as for four of the LRM 
conspiratorial cases discussed above, assuming massless leptons in the final state. As can be 
easily seen here the predictions of the four LRM cases are quite distinct from those of the 
SM even though they all yield the same prediction for B(b —>■ 37). 




Figure 3: Differential decay distribution and lepton forward-backward asymmetry for the 
decay b — > s£ + £~ in the SM(solid) and four representative models in the LRM parameter 
space which yield the SM value for the b —>■ 57 branching fraction and satisfy all other existing 
experimental constraints: Vl = Vr (square- dot), A(l)(dash), A(3)(dot) and D (2) (dash-dot). 



Other observables, such as the longitudinal polarization asymmetry of the r's in 
b — > st + t~ , can be obtained in a straightforward fashion from the expressions above and an 
L — > R augmentation of the expressions provided by Hewett|4]]. These are the only b — > s£ + £~ 
observables we will make use of in the analysis below; the possibility of employing |13[] the 
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transverse polarization of the r's has been neglected. In a similar spirit we ignore the possible 
information that one could gain from future photon polarization measurements in b — > 57. 

4 Analysis Results 

The essential aspects of our procedure can be found in the work of HewettQ] and we closely 
follow the discussion of this author. The analysis makes use of the following observables. 
For the process b — > s£ + £~ , we consider the lepton pair invariant mass distribution described 
by dB/ds, and the lepton pair forward-backward asymmetry, A(s), for £ — e, fi, and r, as 
well as the tau longitudinal polarization asymmetry P T (s). We will neglect the \x mass for 
simplicity and directly combine the e and fi samples. We also include B(b —>■ 37) to this 
list of observables. The lepton pair invariant mass spectrum is divided into 9 bins which are 
distributed as follows: 6 bins of equal size, As = 0.05, are taken in the low dilepton mass 
region below the J/ip resonance, 0.02 < s < 0.32, and 3 bins are in the high dilepton mass 
region above the if)' pole and are taken to be 0.6 < s < 0.7, 0.7 < s < 0.8, and 0.8 < s < 1.0. 
By using this set of bins we completely avoid the regimes where both long distance and 
resonance contributions are clearly important. 

Our analysis proceeds as follows. For a given conspiratorial choice of Vr derived 
above we use Monte Carlo techniques to generate binned "data" associated with the above 
quantities. For the b — > s£ + £~ observables we assume that the errors will remain statistically 
dominated, while for b —>■ 37 we assume a purely systematically dominated error of 7% 
arising from both experimental and theoretical uncertainties. These distributions are then 
generated for an integrated luminosity of either 5 x 10 7 or 5 x 10 8 BB pairs; these correspond 
to the expected total luminosity of a couple of years of running at future I?-factories on the 
T(4S) and at the LHC, respectively. For b — > s£ + £~ the number of events per bin is given 
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by 



bin 



-"I 

J s 



Smax 

ds 



ds , 



(6) 



and the integrated average value of the asymmetries for each bin is then 

{A) h ^^-— / A(s) — ds, 

-''bin -' s mm 

where dB /ds can be obtained from the double differential expression above. 



(7) 
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Figure 4: Comparison of a typical Monte Carlo event sample generated with the right-handed 
mixing matrix -8(3) for (a) the number of events per bin and (b) the bin averaged asymmetry, 
e and fi final states are summed and the histogram represents the best fit to the data. A 
sample of 5 x 10 s BB pairs has been assumed. 



Once the data is generated for each model we then perform a three dimensional x 2 
fit assuming only the usual three coefficients are present, i.e., Ctl,9L,wl(^)- This is done 
according to the usual prescription 



E 

bins 



5Q t 



(8) 
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where Q° hs , Qi(Cl), 8Qi represent the "data", the result of the expectations for a given 
set of Ctl,9L,iol{ij) values, and the error for each observable quantity Qi. The Ctl,9L,iol 
are then varied until a y 2 minimum is obtained. Note that there are 27 bins of data in 
addition to B(b — > 57); allowing for the three free parameters CVl,9L,iol means that the fit 
has 25=27+1-3 degrees of freedom. If the SM were realized, we would expect such a fit 
to yield the SM values of C 7 l,9L,iol, within errors, with a good x 2 — 25 as was shown by 
Hewett[|]]. If new physics is present but the operator basis is not extended we would again 
expect a comparably good x 2 fit at values of Cjl$l,iol which would now exclude the SM at 
some confidence level. In our case with an extended operator basis we would hope that the 
best fit obtained by varying these three coefficients alone is not very good thus demonstrating 
that the three parameter fit is insufficient. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 



2/2 2/2 

s=q /m b s=q /m b 

Figure 5: Same as the previous figure but now for matrix A(l). 

We now give a few examples of this type of analysis. We first consider a typical 
Monte Carlo data sample generated for the matrix -8(3), assuming 5 x 10 8 BB pairs are 
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produced, and examine the results of the best fit. This is shown in Fig. ^ for the two 
observables N and A for combined e and [i data samples. Note that in the low s bins 
the best fit underestimates (overestimates) the number of events at the low (high) energy 
end and generally underestimates the asymmetry. Combining all of the observables this 
particular "data" set leads to a x 2 /d.o.f of 187.1/25. We note that for 25 d.o.f., x 2 values 
of 37.65(44.31,52.61,60.14) correspond to probabilities of a consistent fit of 5(1,0.1,0.01)%, 
respectively. Thus it is clear that for this particular sample the "data" are not consistent 
with the assumption of only three active operators and that an extended operator basis is 
required. 

Figures |^ and ^| show examples of generated data and the corresponding best three 
parameter fits for the right-handed matrices A(3) and D(2). In both these cases the best 
fits are simply incapable of obtaining the correct shape presented by the data. For these 
two cases the y 2 /d.o.f. from these fits are found to be 1187.2/25 and 764.1/25, respectively, 
for a sample of 5 x 10 8 BB pairs. Thus we again see that these particular data sets are not 
consistent with the three Wilson coefficient fit hypothesis. 

Of course, to really ascertain if the new physics of the extended operator basis should 
be visible we need to generate many sets of data for each of the Vr assumptions and determine 
the fraction of the time that the resulting x 2 values exceed those listed above for each of the 
fixed probabilities. To be specific we generate 1000 sets of data for each of the models above 
and perform the x 2 fitting procedure for each data set. For very large data samples, i.e., for 
5 x 10 8 BB pairs, we find that for all of the above choices of Vr we obtain fit probabilities 
below 0.01% almost 100% of the time. This means that it would be quite clear in this case 
that an extended operator basis is required. For smaller data samples, i.e., for 5 x 10 7 BB 
pairs, the results are much more model dependent and are given in Table 2. It seems that 
Vl = Vr represents the worst case scenario. (To test the stability of our results we generated 
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10000 sets of data for this case and confirmed the results shown in the Table.) Obviously it 
is quite important to decide at what level of probability one is willing to exclude the three 
operator fit before drawing any conclusions about an extended operator basis. However, it 
is clear that even for very low probabilities, ~ 0.01%, we see that the three parameter fit 
will fail on average a reasonable fraction, ~ 45%, of the time. 
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Figure 6: Same as the previous figure but now for matrix D(2). 



5 Discussion and Conclusions 

The inclusive rare decays b — > S7 and b — > s£ + £~ have been and will continue to be subjects 
of intense interest since they offer unique opportunities to probe for new physics beyond 
the Standard Model. Both these decay modes are quite clean theoretically and future B- 
factories will produce large event samples in both cases, allowing for in-depth studies of 
their associated observables. As we know, in the SM and in many of its extensions these 
decay observables can be expressed in terms of only three a priori unknown parameters, 
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Matrix P = 5% P = 1% P = 0.1% P = 0.01% 



v L = v R 


~ 39% 


~ 23% 


~ 11% 


~ 6% 


A(l) 


~ 100% 


~ 100% 


~ 100% 


~ 100% 


A(3) 


~ 90% 


~ 76% 


-51% 


~ 30% 


5(2) 


-47% 


~ 28% 


~ 15% 


~ 8% 


5(3) 


~ 92% 


~ 80% 


~ 56% 


~ 36% 


D(2) 


~ 100% 


~ 100% 


~ 100% 


~ 99% 



Table 2: Fraction of the time the effects of the extended operator basis are observable, with 
different choices of the consistency of fit probabilities, for each of the Vr choices described 
in the text assuming a sample of 5 x 10 7 BB pairs. These results are based on 1000 Monte 
Carlo data samples for each of the models. 

corresponding to the values of the Wilson coefficients of three operators at the scale \x ~ m b . 
The global fit approach provides the best model independent technique for obtaining the 
values of these coefficients at the low scale which can then be compared with the expectations 
of a given model. In this paper we have demonstrated, using the Left-Right Symmetric Model 
as an example, that with the statistics available at future S-factories it will be possible to 
observe the rather unique situation where this global fit to the canonical three coefficients 
fails. This result would tell us that not only does new physics exist beyond the SM but that 
this new physics requires an extended operator basis. 
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